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The nucleon magnetic moment in the e-regime of HBChPT 



Joseph Wasent 
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The nucleon magnetic moment is calculated in the e-regime of Heavy Baryon Chiral Perturbation 
Theory to order e^, using the method of collective variables to integrate nonperturbative pion zero 
modes. Contributions containing multiple sources of zero modes enter, allowing for charge-carrying 
zero mode pion fields that connect the sources. The result of this calculation will allow for lattice 
QCD calculations involving nucleons to systematically extract the leading low energy coefficients of 
Heavy Baryon Chiral Perturbation Theory with electromagnetic interactions. 



I. INTRODUCTION 

As the only known way of computing nuclear observables directly from QCD, lattice QCD has received steadily 
increasing attention in recent years. However, limited computational resources have required calculations to use 
unphysically large quark masses and severely restricted computational volumes. Some recent calculations have used 
I ' pion masses that are still a factor of two greater than the physical value, with a spatial lattice volume no greater than 
D _ a few Fermis on a side. Both of these restrictions necessitate the use of Chiral Perturbation Theory (ChPT) in finite 



volume to extrapolate the large mass, finite volume lattice results to physical mass, infinite volume results that can 
be directly compared to experiment. 

is used when calculating baryon properties. When 
electromagnetism is included, simple extensions to the theory that maintain the chiral symmetry are needed and are 

m 



" ^ ■ be directly compared to experiment. 

■ Heavy Baryon Chiral Perturbation Theory (HBChPT) jl, 2| i 



7-H given in refs. [3|, Ij, |5|, |6[ . Typically the small expansion parameters in HBChPT (in what is known as the p- regime) are 

o 

p/niB, vl s-iid m7r/A^ where "p is the typical momentum and is the chiral symmetry breaking scale, typically 



of order 1 GeV[7]. However, as the zero momentum pion propagator is given by Xjm^^V (where V is the spacetime 
volume), for small quark masses the pion zero momentum modes will become enhanced relative to nonzero modes and 
become nonperturbative near the chiral limit ^ . The regime where this occurs is known as the e-regime, and in this 
regime a different counting scheme is used which accounts for the nonperturbative zero modes [T 



In the e-regime the counting parameter is defined by e ^ 27r/A^L ~ 27r/A^/3 and e^ ~ m-^jh.^ where L and /3 
are, respectively, the spatial and temporal dimensions of the box one is performing the calculation in and m„ is the 
pion mass. The inclusion of decuplet baryons in the calculation further requires the definition ^ A/A^, with 
A the decuplet mass splitting. The nonperturbative zero modes are integrated over using the method of collective 
variables ^ and baryons are included in the collective variable framework using the method from ref. jisl . 
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The magnetic moment has been the subject of intense study by the lattice QCD community, including theoretical 



studies in finite [igI. [tI [ist and infinite volume 

29|, |30|. Most of these have involved the p-regime. With faster computers becoming available, lattice theorists will be 
faced with a choice: simultaneously decrease quark masses and increase lattice sizes so as to remain in the p-regime, 
or keep lattice sizes constant and push to lower and more physical quark masses leading directly into the e-regime of 
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221 as well as lattice calculations 
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ChPT. To facilitate the analysis of data from lattice calculations in the e-regime this work focuses on the calculation 
of the magnetic moment of the nucleons in HBChI 
zero momentum modes using the method from ref. 



of the magnetic moment of the nucleons in HBChPT in the e-regime to order e^, with exact integration of the pion 

f. Q- 



II. HEAVY BARYON LAGRANGIAN WITH ELECTROMAGNETISM 



The low energy HBChPT Lagrangian that is consistent with spontaneously broken SU{'2,)l® SU{2)ii is, at leading 
order 0,0: 



/■2 f2 

Co ^ Niv ■ VN - T^iv ■ PT^ + AT^T^ + —Ti'ld^.T,'' d^^T.] + A— Tr[m,I]t + h.c] 

8 4 

+2ff^7V5^^^iV + g AN[f ''''''' Ai,me,d + h.c] + 2gAAnS^A^T-' 
with the nuclcon fields N, the Rarita-Schwinger fields T^^ describing the A-resonances, and the definitions: 



(1) 



= exp 



/ 



M 



7r7V2 7r+ 
n- -7r7V2, 



2?^ = d^ + V^. 



(2) 



where the constant / = 132 MeV, A is the decuplet mass splitting, and the matrix rriq is the quark mass matrix. The 
pion fields are encapsulated in the matrix M. The Rarita-Schwinger fields contain the A-rcsonances according to: 



rpiii ^ ^++ rpii2 ^ -i=A+ r^2^ = ^A° T^^^ = A" 



while the nucleons are an SU(2) vector given by 



N : 



(3) 



(4) 



The three couplings given in eqn. ([T]) {g\, gAN, and (?aa) are the infinite volume, chiral limit couplings between 
baryons and pions. is the covariant spin vector and is the heavy baryon four velocity with — 1 (typically 
calculations arc done in the baryon rest frame with = (1,0)). 



To calculate the magnetic moment in HBChPT, nucleon and decuplet magnetic moment operators are necessary 
along with electromagnetic couplings. Specifically, Feynman graphs that have a component proportional to either 
NF'^'^a^jyN or NF'^'^a^,yT^N are sought. The lowest dimension operators that achieve this are given by the Colenian- 
Glashow SU(3) relations (adapted to two flavor HBChPT) 



LSion opera 

Cmagi = ^ F^, {^idTt {ija^''' ^b) + ^fTt {^a^'' rl^^b)) 



4m 



F^, {^IoN<J^'''N + ^IlN<J^'•'T|+N) (5) 



where the V' are octet baryon spinors and r^_|_ = ^CC^^^C + C'^'^C^)- The decuplet and nucleon-decuplet transition 
magnetic moment operators are determined by finding SU(3) invariants as in ref. 6| and adapting these to the two 
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fiavor theory. These contributions give the additional terms in the Lagrangian 

TUB zniB 



Electromagnetism is included into HBChPT by substituting in eqn. ([T]) for V^, A^, and c)^S according to 

V V^^ + ^ieA'^i^^Q^ + ^Q^^) 

^ A^- ^leA^ (CQ^t _ ^tg^) 

— > P'^E = a^E + zey^'^ [Q,I]] (7) 



and in eqn. ((6]) for Q with: 



Q ^m^ + ^^QO- (8) 



Higher order Lagrangian terms that will also be important include: 



IvfiB '2,'mB 



The chirally covariant derivatives V contain y terms which will not be important at the order considered and so 
are dropped. The coefficients of these terms can be determined by reparametrization invariance 31| or by matching 



to the relativistic theory in the path integral '32|. The resulting terms in eqn. ^ can be included in the baryon and 
decuplet propagators. Following an expansion in I/tub one will then recover the terms equivalent to an insertion of 
this kinetic energy operator. Other l/m^ operators exist beyond these kinetic terms, but they will not enter until 
order e"* or higher, and so are not important for this calculation. 

III. ZERO MODE INTEGRATION AND RESULTS 



In calculating the nucleon magnetic moment in the e-regime the pion zero modes will be determined nonperturba- 
tively using the method described in ref. [15]. Many of the graphs contributing to the magnetic moment will have a 



zero mode structure similar to that of the graphs in ref . [15l | , thus simphfying the calculation. In making the change 
of variables to allow the zero mode integration, a Jacobian factor arises which adds another term (of order e^) to the 
Lagrangian|9|: 

'1 



C-JacoUan = ^pp]^3 [^'^0 + '^+^- ] (10) 



where the n fields are nonzero mode pion fields (following the conventions of ref. [15|). This has the form of a shift 
in the squared pion mass of Am^ = 8/3/^/3L^. When applied to the pion propagators in the order e graphs the 
contribution comes in at order e^, however because it has the form of a mass correction it will be placed in that 
position. 

The leading and next to leading order contributions to the nucleon magnetic moment are given in fig. ([T]). The 
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(a) (b) (c) 

FIG. 1: (a)Leading and (b,c)next to leading order contributions to the nucleon magnetic moment. 

contribution from fig. ((T^) stems from two different Lagrangian terms in eqn. ([6]), one proportional to /iq and one 
proportional to /ii. As the term proportional to /io does not contain any pion field operators it cannot have any zero 
mode structure, while the term proportional to fii has zero mode structure similar to an insertion of the axial current 
in ref. [fsj l . The two order e graphs will have zero mode contributions from the electromagnetic vertex only. At order 

the graphs in fig. ([2]) are found to contribute, where figs. ((2^) and (l2t) are the leading wavefunction renormalization 
contributions. Again these contributions are very similar to those found in the calculation of the axial charge, and as 
such have the same zero mode structure. 

At order most of the new contributing diagrams are formed by attaching a simple pion loop at every vertex in 
an order e diagram, as in fig. ([3]a,-d). In some cases (figs. ^) and ([3ji)) this does not create any new zero mode 
structure, as all of the zero mode information comes from the electromagnetic vertex. However, in fig. ([3^) and ^p) 
the addition of the pion loop to the nucleon-pion vertex creates a new source of zero modes in addition to those 
from the electromagnetic vertex. The zero mode functions from contributions of this form separate into two groups 
of diagrams. The first group (characterized by the function A(s) as defined in the appendix) contains intermediate 
nonzero mode pious where charge is conserved separately within both the zero mode and nonzero mode sectors of the 
theory. The second group (characterized by the function B{s)) contains intermediate zero mode pious which carry 
electric charge from one vertex to the other. The last two diagrams in fig. ([3]) are not as simply constructed, but are 
straightforward to compute. Figure ([3]3) is the lowest order two loop construction (the electromagnetic vertex can go 
on either pion loop with an equal contribution) while fig. ([3f ) is an insertion of the pion kinetic energy operator into 
an order e diagram. 
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(g) 



FIG. 2: Contributions to the magnetic moment at order e^. Graphs (e) and (f) are wavefunction renormahzation diagrams. 




(e) (f) 
FIG. 3: Examples of contributions to the magnetic moment at order e^. 



Once the various zero mode functions have been determined, the remaining finite volume nonzero mode Feynman 
diagrams can be written in the standard way as a four component sum over integers with the zero vector specifically 
removed. Using the Abel-Plana formula, the fourth component of this vector sum can be written as an infinite volume 
integral plus finite volume corrections. As shown in ref. due to the nature of the HBChPT propagators, these 
correction terms will be negligible until at least order and so can be ignored. Using this fact, the diagrams can be 
determined by integrating over the time component of the momentum and then summing over the remaining spatial 
components. To perform the sums that arise, one can explicitly expand them in powers of e to the order desired, 
resulting in purely numerical sums (expressed below as coefficients defined in the appendix). 
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Calculating the contribution from each of the graphs given in figs. ((T][3]) results in: 

jl = A0 + A1+A2 + A3 (11) 

Ao = Mo + MiT^Q(s) (12) 

^mB{g\f C2 C4i(TO^(s) + 8/3/2/3L3) 1 3ci , 1 co^i , IQitibqIn 3^, ^ ^ C2 



B{9Ar ^ / C2 C4-L(TO^(,SJ /jj.-J 1 .^Ci i Cp \ IbjTlBff^jy 3 / C2 

3/2 ^^"^^^vS^'i 167r4 ^ 2tob 8^i2 + 4^2 ^3 y + 27/2 ^^ ^8^22, 



3Ac: 



■3 



(A2-(m2(g) + 8/3/2/3L3))£c4 ^ _ Ac^\ 1 cp 



327r3 167r4 2tob VSTri^ 2Tr^L I 4m% 



(13) 



= Q(^)^-^(^) + 372^ - 



27 /2 '"^^ V47ri2 87r2L 2mB J 27r/2L2 

- (a^o + Mir Q(,s)) + ^—^ J + 3 7^ (^i^ ~ 4^ + 2;;^^ J J (^^^ 



8/4713X3^^^-4^ 81 9^Aivv..Ayy -^v y 97r3/4L3 Sln^f^L^ 

C3™B 3^/ \ \4 , 16/ p n2 2 84 56 q 2 400 2 2 \ /^r\ 

rr^Qis) {9a) + ^(5a) " ^5ajv " ^ffAffAA^AJV - ^SaatSaa (15) 



367r3/4L3' ^^"^ ' 3 V^A; ^AiV g-^^AiV g^^A^^-^AiV ^29^ 



with 



^ A' ^ N , 8 ^, / A~v/A2-(m2(.)+8/3/2/3L3) >^ ^ m2 (.) + 8/3/2/3L3 ^ 

where the subscript refers to the order in e of the contribution. The zero mode functions Q{s), T{s), A(s), and 
B{s) as well as the coefficients cp, ci, C2, C3, C4, and c^i are defined in the appendix. The pion mass also has a zero 
mode contribution that has been previously calculated 8| which effects the replacement m'^ —^ ■m^{s). The zero mode 
structure of this replacement is the same as that from the tadpole diagram in fig. ([2^) and indeed m'^{s) = ra^Ti^s). 
Note that the infinite volume portions of the order e contributions enter at order , and are the first two terms in /i2. 

IV. DISCUSSION AND CONCLUSION 

We have computed the nucleon magnetic moment in the e regime to order e^. In doing so, two new zero mode 
functions have been calculated, A{s) and B{s). Both of these functions have intriguing features, as each stem from 
diagrams with two separate vertices that produce zero modes. The zero mode effects contained in the function A{s) 
are generated from zero mode fields that connect the vertices but do not carry charge, as shown in fig. (U^). This is 



in contrast to the zero modes encapsulated in the function B(s), as they connect and carry charge between the two 
vertices, as shown in fig. (|Dd). When writing down contributing diagrams where the zero modes have been removed 
and placed in the function B{s) the diagrams will appear to violate local charge conservation. In fact, charge is 
locally conserved when the zero modes are included. The form of the function B{s), as shown in fig. (O, is easy 
to understand. As s — > oo the zero mode fluctuations become damped and contributions of this type disappear due 
to charge conservation in the nonzero mode pions, while as s — > the charged zero modes will become increasingly 
important. 





(a) (b) 

FIG. 4: Examples of the type of zero mode pions that contribute to the functions (a)yl(s) with charge neutral zero mode 
configurations (grey blobs) connecting each vertex and (b)_B(s) with charge carrying zero mode configurations connecting each 
vertex. 

The behavior of the zero mode function A{s) is unique because it changes sign as s ^ 0, raising an interesting 
possibility. For a given choice of spatial volume and pion mass, the temporal dimension of a lattice calculation could 
be chosen such that A(s) was zero, thus removing these contributions from the calculation. In addition, previous 



studies (see refs. 
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33|,|3J]) have discussed the changes that occur in the coefficients ci and C2 in asymmetric spatial 
volumes. Specifically, these coefficients can be either enhanced or made to vanish with the correct choice of spatial 
dimensions. If one was to combine these spatial dimension effects on the coefficients with the temporal dimension 
effects on the zero mode functions one could design a lattice with a significant degree of control over the contribution 
of certain graphs. 

Due to the increased complexity of lattice calculations involving electromagnctism over lattice calculations with 
pure QCD, the lattice theorist may be tempted to forgo working on quantities such as the magnetic moment in favor 
of quantities such as the axial charge. However, the electromagnetic properties of the nucleons are some of the most 
precisely measured quantities in physics, and as such an important test of lattice QCD will be the reproduction of these 
values. Because of this, lattice calculations of the nuclcon magnetic moment in the e-regime will be of importance in 
the coming years. Furthermore, configurations used to calculate quantities such as the axial charge can also be used 
to calculate the magnetic moment. In fact, calculating both quantities simultaneously and performing the appropriate 
correlated statistical analysis on each would be a powerful tool for determining the infinite volume low energy coupling 
coefficients g^, gAN, and ^aa sls well as the magnetic couplings /ioj A^ii A'Ci ^^nd /it- With the above results for the 
magnetic moment and previous results [fsjl for the axial charge, both quantities can be coherently analyzed in the 
e-regime. 
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APPENDIX A: ZERO MODE FUNCTIONS 



The zero mode functions (with s = ^m^f V) calculated for the contributing diagrams are: 



Tis) 
A{s) 
Bis) 



3 I + /i(2.) 
/o(25) 1 
/i(2s) s 

11 1 Jo(2,s) 2s/o(2s) 

6 3s * 3/i(2s) h{2s) 



3 J_ 

5 ^ 5l 



8s , /o(2s) ^ 8s/o(2s) 
"5" 



(Al) 
(A2) 
(A3) 
(A4) 



5/1 (2s) 5/1 (2s) ■ 

Of these functions Q{s) and T{s) have been calculated previously in refs. PjIi^. The functions A{s) and B{s) are 
new to this paper and are generated from contributions of the type shown in figs. ((3^) and ^p). These functions are 
plotted in fig. 




0.2 



FIG. 5: The form of the four zero mode functions used in the calculations. 
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APPENDIX B: SUMS 

For the calculation of the dia gram s given above several sums will be important. Many of these sums have been 
calculated previously 
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37[ and the divergent sums are defined through dimensional regularization: 
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= -2.8372974 
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fi#0 







CO = 2.1 = -1- (Bl) 

The sum that needs to be evaluated for the two loop diagram is more complicated, but can be evaluated using the 
same general techniques: 

4 = E 1^ 1^1,1^1^ 1^1,, + T^TTTr^r-^^] = 2.722962. (B2) 
^ ^ \ 2 |n|(|m| + |n|2(|m| + |n|)2 / ^ ' 
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